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Abstract 

The quantum Hall effect in the three-dimensional anisotropic tight-binding 

electrons is investigated in the field-induced spin density wave phases with a 

magnetic field tilted to any direction. The Hall conductivity, axy and axz, 

are shown to be quantized as a function of the wave vector of FISDW, while 

Gyz stays zero, where x is the most conducting direction and y and z are 

perpendicular to x. 
73.40.Hm, 75.30.Fv 
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I. INTRODUCTION 



Quantum Hall effect has been observed in organic conductors, (TMTSF)2PF6 and 
(TMTSF)2C104§H1 in the field-induced spin density wave (FISDW) phases. These organic 
conductors are quasi-one- dimensional conductors and described by a tight binding Hamil- 
tonian with anisotropic hopping matrix elements, t;,/ta ~ 0.1 and tcAa ~ 0.003. The Fermi 
surface consists of two warped planes in the absence of a magnetic field. If the nesting of 
the Fermi surface is perfect, the staggered susceptibility corresponding to the nesting vector 
which connects two planes of the Fermi surface diverges as temperature becomes low. In 
such case the spin density wave (SDW) is stabilized by the repulsive interaction of electrons. 
The nesting of the Fermi surface is imperfect in general due to transverse hoppings. For 
the imperfect nesting case the transition temperature of SDW depends on the imperfectness 
and the strength of the electron interaction. Since the transverse hoppings are enhanced 
by pressure, the SDW state is suppressed as the pressure is increased. The divergence of 
the staggered susceptibility recovers when the magnetic field is applied@, and the successive 
transitions to FISDWs occur as the magnetic field is increasedi^. One can explain the 
stability of FISDW by the fact that the energy is lowered when a gap at the Fermi surface 
is opened due to the combined effect of the electron interactions and magnetic field. 

A magnetic field required to be one flux quantum, 0o = hco/e, per electron in each plane 
for these organic conductors is extremely large (~ 10^ Tesla). In other words the quantum 
Hall effect in organic conductors is observed in much smaller field than that is necessary 
to fill electrons only in the lowest few Landau levels. Moreover, the Hall conductivity is 
observed to change its sign in some region of the magnetic field0. The Quantum Hall effect 
in the presence of FISDW has been studied theoretically0Hll. The quantization of the Hall 
conductivity is explained by using the Streda formulaS or the general theory of the Hall 
effect in the periodic syst These authors studied the case that the magnetic field is 

perpendicular to the conducting plane. 

Due to the hoppings along the third direction new phenomena are expected when the 
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mametic field is tiltecfl. In three-dimensional lattice there exist three fluxes, 0a, 06 and 0c 
per unit area in each plane, while only one flux per unit cell is present in two-dimension. 
Indeed, many cusps of the magnetoresistance have been observed as a function of a direction 
of the mae; netic fleld0'0. FISDW in three-dimension in the tilted magnetic fleld or in non- 
orthogonal lattice has been studied^. The Hall conductivity a^y and predicted 
to be quantized. These values are argued to be fractionally quantized@'0. These authors, 
however, studied the case that the magnetic fleld is perpendicular to one of the axis, say a, 
when two fluxes 0b, and 0c exist instead of three in the general direction of the fleld. 

On the other hand the general formula for the quantum Hall effect in the periodic system 
in three-dimension was given by Halperin^l and Kohmoto et a/.ii. They have shown that 
the conductivity tensor is quantized when the Fermi energy is in the energy gap. 

In this paper we study the Hall conductivity in the anisotropic three-dimensional lattice 
with the magnetic fleld tilted to any direction. In section the conductivity tensor is 
calculated for the non-interacting electrons in the orthorhombic lattice. The quantized value 
is given in the perturbation in t^/ta and tc/^a- The quantum Hall effect in the presence of 
FISDW is given in section |ITT[ The generalization to the triclinic lattice is done in section 

II. NON-INTERACTING CASE 

In this section we study the anisotropic tight-binding electrons on the orthorhombic 
lattice in the magnetic fleld, 

{i,j)a,cr {i,j)b,o- (i,j)c,(J 

where cj,^ and q^o- are creation and annihilation operators of electron with spin a at site 
ta, th and tc (ta ^ U, tc) are the hopping matrix elements along a, b and c directions, 
respectively, and 

27r P 

(j)ij = — A dl, (2.2) 
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where A is a vector potential. In this section and section |T| the a, b and c axes are assumed 
to be orthogonal each other and parallel to the x, y and z direction, respectively. 
The uniform magnetic field, B, is applied in any direction, 

B={B,,By,B,). (2.3) 

The fluxes through the unit area perpendicular to the a, b and c axes are given by 

0^ = bcB^, (2.4a) 

(py = caBy, (2.4b) 

0. = abB„ (2.4c) 

where a, b, and c are the lattice constants. 

First we consider a rational flux case, i.e., the flux through a unit area in each plane is a 
rational number, {(j)x/4>o,4>y/4>o,4>z/4>o) = {Px/qx,Py/qy,Pz/qz) with mutually prime integers 
Pa and Qa {a = x,y or z). Then the Hamiltonian in the presence of the magnetic fleld 
is described as a generalization of the problem in two dimension studied by AzbelS and 
Hofstadterii to three dimension0'&^. The energy spectrum for the rational flux case can 
be obtained numerically as eigenvalue of Harper equations. The size of the Harper equations 
is the least common multiple of qx, qy and g^, which we deflne Q. The volume of the magnetic 
Brillouin zone is 1/Q of the Brillouin zone in the absence of a magnetic fleld. For a fixed 
momentum k in the magnetic Brillouin zone there exist Q eigenvalues and Q eigenstates. 
By varying the momentum we get Q bands. In two dimension these bands do not overlap. 
Although the bands may overlap in three dimension, it has been shown that the energy gaps 
due to a magnetic fleld exist in the wide range of parameters00. We can use the extended 
zone by unfolding the magnetic Brillouin zone to the Brillouin zone in the absence of a fleld. 
Then the energy is obtained uniquely for each momentum. 

The quantum Hall effect is expected in the three-dimensional case when some of the 
energy bands are completely fllled and the other bands are emptyi'S'ili. The conductivity 
tensor for a fllled band is given as^iiii 



where the factor of 2 comes from the spin degrees of freedom, a and jS are x, y, or z, > is 
the wave function for the filled band and the integral is performed in the magnetic Brillouin 
zone. Kohmoto, Halperin and WuS have shown that if the Fermi energy lies in the energy 
gap, the conductivity tensor is described as 

Ca/S = 2;- — -realS^Gy, (2-6) 

znn 

with the vector in the reciprocal lattice, 

G = —{la'Ga' + h'Gb' + Ic'Gc') (2.7) 

where Ga', Gb' and Gc' are the fundamental reciprocal vectors satisfying B-Ga' = B-Gf,' = 
and la', lb' and Ic' are integers. These integers are the first Chern numbers on the tori obtained 
by slicing the three-torus of the magnetic Brillouin zone. Since we can take any fundamental 
reciprocal vectors, G is written as 

G = -ilaGa + lbGb + lcG,), (2.8) 

where Ga = (27r/a)x/|x|, G;, = (27r/6)y/|y|, Gc = (27r/c)z/|z| for the orthorhombic lattice 
and la, lb and Ic are integers. If the direction of the magnetic field is changed infinitesimally, 
G stays constant as long as the Fermi energy lies in the energy gap. Although Ga', Gb' and 
Gc' depend on the direction of the magnetic field, G^, G;, and Gc do not depend on the 
magnetic field. As a result, while la', h' and Ic' are not constant, la, h and Ic are constant 
for the infinitesimal change of the magnetic field. 

If the Fermi energy lies between the rth and r + 1th band from the bottom of the energy, 
the conductivity tensor is given as the summation of the contributions from filled n bands. 
The contributions from the filled bands cancel each other except that from the upper energy 
gap of the rth band0. 

In order to calculate the Hall conductivity explicitly, we take the vector potential A as 



A = (0, (B.x - B,z), -^{B,x - B^z)). (2.9) 

Bz 

In the above we have assumed B^ without loss of generahty, since we may exchange y 
and z in the case of Bz = 0. With this vector potential the non-interacting Hamiltonian is 
written 



"^0 = -iaX]2cos(a/c3;)4(k)C(,(k) 

k,(T 

-tbEe"*'''4(k-u)c,(k)+/i.c. 



k,(T 

._ (px<Py 



-t^e"^ E e-*"^^4(k + u')c^(k) + h.c, (2.10) 

k;,(T 



where /i.e. means Hermitian conjugate, 



and 



2Ti(j)z „ 27r0^\ 

u= — 1-,0, — , 2.11 

V a 00 c 00 / 



u' = ^u. (2.12) 



We define the state as 



^,(k) = 4(k)|0>, (2.13) 

with vacuum |0 >. The Hamiltonian mixes ipai}^) with ^/'o-(k±u) and ■?/'cr(k± u'). If (j)x = 0, 
ky and /c^ are constants of motion. If (p^ 7^ 0, kz is not a constant of motion but ky and 



{(f)zkz + 4>xkx) I Y^0^ + 02 are constants of motion. 

The Fermi surface is given by kx = -^kp in the zeroth order perturbation in tb/ta and 
ic/ia- When the condition 

2kF — mzUx-\- myu'^ + s{2t: /a) (2-14) 

is satisfied with integers my, rUz and s, the degenerate states at k^ — kp and k^ — —kp are 
mixed resulting the energy gap at the Fermi energy in the Im^jth perturbation in (tb/ta) 
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and the |my|th perturbation in (tc/ta). The magnitude of the energy gap is of the order of 
|ro(k)|, where 

^"("^-'-fc) it) ' 

= tj£^) (h\ (T^-^l-^ ^K-m.bky+myck.) _ (2.15) 



^^a/ Via/ 

Since the flux is assumed to be rational, u and u' are written as the integer multiple of the 
vector uq as 

u = PzUo, (2.16a) 
u' = PyUQ, (2.16b) 



where integers Py and are defined by 

Py_^(i)y_ 

Pz (l)z' 

The Harper equations are obtained by using the one dimensional basis states 



(2.17) 



^j = ^{k + juo), (2.18) 

with integer j. The number of the Harper equations obtained by this basis set is not always 
Q but the integer multiple of Q, because ipQ is not necessarily the same as i/jq. Thus the basis 
tpj may be overcomplete, although the eigenenergy is obtained correctly by the overcomplete 
basis. In order to get the complete basis for the rational flux case a careful choice of the 
vector potential is necessaryil. The overcomplete basis set corresponds to the larger unit 
cell in the real space or the folded magnetic Brillouin zone in the momentum space. The 
Hall conductivity is, however, obtained even if the integration in Eq. ( |2.5| ) is done in the 
folded magnetic Brillouin zone, since the degeneracy of the bands due to the overcomplete 
bands compensates the folded region of the integration. 

The magnetic Brillouin zone has the volume {27t)'^ / (Qabc) and 2kFaQ/{2n) = r is the 
number of filled bands below the Fermi energy. Then we get from Eq. ( |2.14] ) that niy and 
TTLz should fulfill the Diophantine equation. 



r = sQ + -j^QrUz + -j-Qniy. 



(2.19) 



The wave function for the rth band gets the phase of ro(k) at ~ k^^, i.e. the phase 
of the wave function cannot be defined globally in the magnetic Brillouin zone. On the other 
hand the phase of the wave function can be defined globally in the torus of ky and k^ for 
fixed kx, which is the slice of three-torus by the plane perpendicular to a fixed k^. Therefore, 
the second term in Eq. (|2.5|) is zero by partial integration with respect to ky or kz if (3 = y 
or P = z. Then the conductivity tensor is written asili 



le' 1 



a. 



xy 



c h 2TTi 



le2 1 

O'xz = 27 



b h 27ii J-tt/c 



w/b 

-Tv/b 

tt/c 



dkylK"^ 



d 


^ > 


dky 


d 




^ > 


dkz 



- < ^1 



d 



dk 

- < '^l — 

k^=k- dkz 



1^ > 



y 



1^ > 



kx — k r 



kx — k r 



ayz = 0, 



(2.20a) 

(2.20b) 
(2.20c) 



where kp{kp) is the momentum smaller (larger) than the Fermi momentum. Since the 
wave function changes only in the phase at the Fermi momentum, the integration with 
respect to ky or k^ gives the winding number of ro(k) in the complex plane around zero 
when ky or kz is moved. Using Eqs. (|2.15|) and ( p.20| ), we obtain 

2e2 



xy 



c h 
2e2 



(2.21a) 
(2.21b) 



The conductivity per plane is quantized as axy = 2{e'^/h)mz and azx = 2{e^ /h)my. 

Unfortunately, the energy gap will be very small if \my \ + |m^| ^ 1, which is the case 
for l^^l < 00- 

For the quantum Hall effect, it is not necessary that all of ipx/'Po-, 4'y/4'o and (f)z/4'o are 
integers. The only required condition is Eq. (|2.14| ) . 



III. 3D QUANTUM HALL EFFECT IN THE PRESENCE OF FISDW 

In this section we take account of the interaction 
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^' = f^E<TC^.T<iC,,i- (3.1) 

i 

This interaction is written as 

n' = -U (^^^ JdkJ dk'J dK c|(k)q(k + K)c{(k' + K)cr(k'). (3.2) 

The order parameter of the spin density is defined in the mean field approximation as 

A(K) = -f/^ J dk' < c|(k' + K)c,(k') > . (3.3) 

The interaction term mixes the states ip^{k) and ipi{k + K). As a resuh an energy gap 
at the Fermi surface is opened by the order parameter when = ±2kp. The y and z 
components of the wave vector should be determined to give the lowest energy, which can 
be calculated numerically with further approximation such as the linearization with respect 
to A;Ji0. Since we focus on the quantization of the Hall conductivity, we do not have to 
calculate the values of Ky and explicitly. The Hall conductivity is zero in this case. 

There are other possibilities for opening gaps in the presence of the magnetic field. The 
order parameter may have many components with respect to the wave vectors as in the 
two-dimensional case&B. In three-dimensional case the order parameter has the form 

A(K)= I^n„n.5(K-K„^,„J, (3.4) 

ny,nz 

with the X component of the wave vector K.ny,n^ given by 

= 2kF + ny — -r + n, — — , (3.5) 



where ny and are integers. The vectors u and u' are given in Eqs. ( |2.11| ) and ( p.l2|) . By 



this wave vector the state at k^ = —kp is mixed with that a.t k;j. = kp + i^yu'^ + n^Ux- The 
state at k^ = kp + Uyu'^ + n^Ux is mixed with the state at kx = kp by |?T.j,|th perturbation 
in thexp{ibky)/ta and jn^jth perturbation in tcexp{ickz)/ta. As a result two states at the 
Fermi energy are mixed by a combined effect of SDW and the magnetic field. 

(See Fig. 1.) The gap function is given by 
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r(k) = J2 Tn^n, exp 



i i-rizhky + Uyckz - tttt,^^^ 

V 0200 , 



(3.6) 



where Tny,n^ is a constant of the order of Dny^nA^b / ta)^'^''\tc/ ta)^'^^^ ■ The wave function near 
the Fermi momentum changes the phase due to this gap function. The Hall conductivity is 
calculated by Eq. (|2.5| ) as in the non-interacting case. 

If only one component of the order parameter is dominant, i.e., 

A(K)^Z}„^„,„,„5(K-K„^„,„,J, (3.7) 

the fc-dependent phase change is exp[i(— n^o^/c?; + nyQck^)] and the Hall conductivity is ob- 
tained by Eqs. ( ^.2(J| ) as 

CTxy = "T^^-^o, (3.8a) 
o-zx = I'f^'^yo^ (3.8b) 

and 

ay, = 0. (3.8c) 

Comparing Eqs. - with Eqs. (^) and (^) we find that 

la = 0, (3.9a) 
lb = riy, (3.9b) 
Ic = n,. (3.9c) 

If the order parameter consists of many components with respect to the wave vectors Kny,n^, 
the Hall conductivity cr^y {cxz) is obtained as a winding number of the gap function r(k) as 
ky {kz) is moved from — 7r/6 to n/h (from — tt/c to vr/c). Note that 0^. appears as the relative 
phase between the components of the order parameter, so the Hall conductivity may depend 
on (f)x implicitly. 

If the magnetic field is changed in magnitude or direction by a small amount, the wave 
vectors and the order parameter will change in order to keep the Fermi energy in the same 
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energy gap. As a result the plateaux of the Hall conductivity is realized without the effect 
of a localization. When the energy gap is closed and opens again at the Fermi energy as 
the magnetic field is changed, the Hall conductivity changes. The sign change of the Hall 
conductivity may happen as the direction or the magnitude of the magnetic field is changed. 

The condition for the rational flux is not necessary to get the Quantum Hall effect. The 
wave vector of the order parameter is adjusted to satisfy Eq. (|3.5|) for any direction and 
amplitude of the magnetic field. If 0y/02 is irrational, Uy and are determined uniquely 
for a given wave vector (Kny,n:,)x- On the other hand Uy and are not uniquely determined 
if (f)y/4'z is rational. 



IV. TRICLINIC LATTICE 

In this section we consider the triclinic lattice with a, b and c not being orthogonal each 
other. The fundamental reciprocal vectors, 

(4.1a) 
(4.1b) 
(4.1c) 





b X c 




= 27T 




V 


Gfe 


c X a 


= 27r 




V 


G, 


a X b 


= 27r 



are not orthogonal, where 

i;=|a-(bxc)|. (4.2) 

We take a || x. Hopping matrix elements t^, 4 and tc (|ta| ^ \tb\, \tc\) are assumed between 
the nearest sites along a, b and c axes. The vector potential is taken as Eq. ( |2.9| ) by using 
the orthogonal x, y and z axes. 

The non-interacting Hamiltonian is 

no= -t„;^2cos(a- k)c^(k)c(k) (4.3) 

k,cr 

- tbC^^" e"'^-^ci{k - w)c,(k) + h.c. 

k,cr 

- t,e''^ Yl e"*"'^4(k + w')c.(k) + h.c, 



where 



Ob = ^{hy - ^h){BX - B^K) (4.4a) 
Oc = ^{cy - |^c,)(S,c, - B^c,) (4.4b) 



and 



'^ = ^(^v- ir^^)iBz^ - B.,z), (4.5a) 

<P0 ^z 

^' = -^-T^^y - |^c.)(i?.x - i?,z). (4.5b) 

Note that w || w'. 

The Fermi surface in the zeroth order in tb/ta and tcjta is the plains given by a-k = ^akp- 
The order parameter is defined by Eq. ( ^.3] ) with the factor v/(2ttY instead of a6c/(27r)^. 
As in the previous section, the energy gap opens at the Fermi momentum by the effects of 
the order parameter and the perturbation in tf,/ta and tc/^a, if the order parameter has the 
form 

A(K)= Y: I^„„n.5(K-K„,,„J, (4.6) 

nt,nc 

where and ric are integers and the wave vector K„j_^„^ satisfies the condition 

^nt.ric ■ ^ = 2akF + rifow' ■ a + n^w ■ a 

o ; , o "CyE^ - ac^By abyB^ - ab^By 

= 2akF + n;,27r — - nc27r — 

00 00 

= 2akF + n^,2^^— + 71^27^—. (4.7) 
00 00 

In the above (ph and 0c are the fiuxes through a unit area in a-c and a-h planes, respectively, 

0c = a ■ (b X B) (4.8a) 
0b = -a ■ (c X B) (4.8b) 

The phase of the wave function is not defined globally as in the previous sections. In the 
perturbation in t^/ta and tc/ta-, the mismatch of the locally defined phase for the wave 
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function of the highest occupied band is attributed to the phase of the gap at the Fermi 
momentum. The wave function changes the phase as the momentum is moved at the Fermi 
momentum as in the previous section. We consider the case that one component is dominant 
in Eq. (|4.6|). Then the k-dependent part of the phase change is 



exp [i {—rich ■ k + nf,c ■ k)] . 



(4.9) 



The conductivity tensor is calculated by 



- < SI 



_d_ 



1^ > 



1 



d 



1^ > 



< ^1 



_d_ 

dk. 



a. 



yz 



0. 



(4.10a) 

(4.10b) 
(4.10c) 



Because of the non-orthogonality of a, b and c, the slice of the magnetic Brillouin zone 
by the plane perpendicular to ky or is not the torus. As a result, a^y and a^x are not 
quantized to be integer but they are quantized as 

o2 



e a 
a^y = 2—[ncby - UbCy)- 



nM 



ribCy 



6 a 
n V 



(4.11) 



(4.12) 



In the above we have used the facts that a || x and the area of the Brillouin zone in 
ky-kz plane is {2iT)^a/v. Comparing these result with the general expression of the Hall 



conductivity in three-dimensional periodic system (Eq. (|2.8| )), we get 



= 



h = rib 



(4.13a) 
(4.13b) 
(4.13c) 



13 



V. CONCLUSION 



When the magnetic field is tilted to any direction, there is a possibility for the new wave 
vector for FISDW. We have shown that the Hall conductivity a^y, Czx and cr^^ are quantized 
as a function of the wave number of FISDW. When the lattice is orthogonal, both axy and 
a^z per plane are quantized as integers and stays to be zero as long astb/ta and tc/ta are 
treated in perturbation. Since the quantization is realized with respect to the fundamental 
reciprocal lattice, axy and azx is not quantize as integers if the a, b and c axes are not 
orthogonal each other. Even in that case the quantum values of axy and a^x are given as a 
function of the wave vector of FISDW, which is characterized by integers. 
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FIGURES 

FIG. 1. Schematic diagram of opening the gap at the Fermi momentum. The right-arrow 

labeled by -Do,o indicates the order parameter of SDW with Kx = 2k p by which two states at 
kx = ifci? are coupled. The right-arrow labeled by Di^i indicates the SDW with wave vector 
Kx = 2kF + Ux + u'x- The left-arrows indicate the coupling of the states caused by the magnetic 
field in the perturbation in tb/ta and tc/ta- 
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